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ABSTRACT 

The energetic, eclipsing millisecond pulsar J1816+4510 was recently discovered in a low-frequency 
radio survey with the Green Bank Telescope. With an orbital period of 8.7 hr and minimum companion 
mass of 0.16 M© it appears to belong to an increasingly important class of pulsars that are ablating 
their low-mass companions. We report the discovery of the 7-ray counterpart to this pulsar, and 
present a likely optical/ultraviolet counterpart as well. Using the radio ephemeris we detect pulsations 
in the unclassified 7-ray source 2FGL J1816. 5+4511, implying an efficiency of ~ 25% in converting 
the pulsar's spin-down luminosity into 7-rays and adding PSR J1816+4510 to the large number of 
millisecond pulsars detected by Fermi. The likely optical/UV counterpart was identified through 
position coincidence (< O'.'l) and unusual colors. Assuming that it is the companion, with R = 
18.27 ±0.03 mag and effective temperature > 15,000 K it would be among the brightest and hottest of 
low-mass pulsar companions, and appears qualitatively different from other eclipsing pulsar systems. 
In particular, current data suggest that it is a factor of two larger than most white dwarfs of its mass, 
but a factor of four smaller than its Roche lobe. We discuss possible reasons for its high temperature 
and odd size, and suggest that it recently underwent a violent episode of mass loss. Regardless of 
origin, its brightness and the relative unimportance of irradiation make it an ideal target for a mass, 
and hence a neutron star mass, determination. 

Subject headings: binaries: eclipsing - ■ gamma rays: stars - - gamma rays: individual 
(2FGL J1816. 5+4511) — pulsars: individual (PSR J1816+4510) — ultraviolet: 
stars 



1. INTRODUCTION 

The 3.2-ms pulsar J1816+4510 was discovered as part 
of the Green Bank North Celestial Cap (GBNCC) sur- 
vey (Stovall et al. 2012, in prep.) — a survey of the sky 
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north of declination +38° at 350 MHz with the 100-m 
Robert C. Byrd Green Bank Telescope — in a pointing 
selected for being coincident with an unclassified Fermi 
7-ray source (a successful search strategy, as shown in 
Hessels et al. 2011, Ransom et al. 2011, Cognard et al. 
2011, Keith et al. 2011, Kerr et al. 2012, and others). 
Shortly after discovery, it was realized that the radio 
data showed evidence for acceleration in an 8.66 hr cir- 
cular orbit with eclipses lasting up to 10% of the orbit 
at 350 MHz (Fig. 1). Eclipsing millisecond pulsars, es- 
pecially those with 7-ray counterparts, are often associ- 
ated with "black-widow" or "redback" systems. These 
systems harbor low-mass companions (< 0.05M© for 
black-widows and ~ 0.2M© for redbacks; Fruchter, Stine- 
bring, & Taylor 1988; D'Amico, Possenti, Manchester, 
Sarkissian, Lyne, & Camilo 2001; Archibald, Stairs, Ran- 
som, Kaspi, Kondratiev, Lorimer, McLaughlin, Boyles, 
Hessels, Lynch, van Leeuwen, Roberts, Jenet, Champion, 
Rosen, Barlow, Dunlap, & Remillard 2009) and have 
been discovered with increasing frequency in recent years 
(see Roberts 2011 for a recent review), often in globular 
clusters. The eclipses are typically long (they can cover 
most of the orbit; Archibald et al. 2009; Hessels et al. 
2011), implying eclipsing regions larger than the Roche 
lobes of the companions, and there are regions of the or- 
bit where the pulsar is seen through ionized plasma that 
delays the pulses compared to the expected ephemeris. 
The basic model for these sources is one in which the 
energetic wind from the pulsar irradiates and ablates the 
companion, leading to long eclipses from ionized material 
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in the systems (e.g., Stappers et al. 2001). The compan- 
ions are usually tidally distorted, filling a significant frac- 
tion of their Roche lobes (Reynolds et al. 2007), which 
along with heating from the pulsar's wind leads to signif- 
icant (often > 3 mag at wavelengths of around 5000 A) 
optical modulation. Such systems are interesting both 
because they provide a probe of the interaction between 
the pulsar's wind and the companion and, ultimately, 
because such systems allow measurement of neutron star 
masses through binary modeling (van Kerkwijk, Breton, 
& Kulkarni 2011). 

Here we report on new and archival optical and ul- 
traviolet data on the counterpart of PSR J1816+4510. 
We use these data, along with archival 7-ray data from 
the Fermi spacecraft, to constrain the nature of the 
PSR J1816+4510 system. The outline of this paper is 
as follows: we first describe the archival optical and ul- 
traviolet data that we used to identify the counterpart to 
PSR J1816+4510 (§ 2.1), and then we discuss new optical 
data from the Wisconsin Indiana Yale NOAO telescope 
(§ 2.2). We fit the optical/ultraviolet spectral energy 
distribution in § 2.3. We then discuss the Fermi 7-ray 
data (§ 3.2) and the Swift X-ray upper limits (§ 3.1). Fi- 
nally, we discuss the implications of our data in § 4, and 
conclude in § 5. 

1.1. System Parameters 

We make use of the radio ephemeris for 
PSR J1816+4510 determined by Stovall et al. (2012): 
position a = 18 h 16 m 35?9314(2) S = +45°10'33'. / 864(2) 
(J2000), binary period Pi, = 8.66 hr, and minimum 
companion mass M c = 0.162Mq (assuming a neutron 
star mass of 1.4M Q , although a somewhat more massive 
neutron star may be likely given likely accretion histo- 
ries; Verbiest et al. 2008; van Kerkwijk et al. 2011). Such 
a companion mass would put it among the "redback" 
class (Roberts 2011), although if the companion is more 
degenerate it might instead contain a He-core white 
dwarf such as that in the PSR J1911— 5958A system 
(Bassa et al. 2006). If the orbit is edge-on, then the 
full semi-major axis is a = 2A6R@ with a Roche-lobe 
radius R L = O.53i? (based on Eggleton 1983). The 
dispersion-measure (DM) distance is 2.4 kpc (Cordes & 
Lazio 2002, for a DM of 38.8 cm" 3 pc), although given 
the high Galactic latitude (b = 24.7°) the uncertainties 
are large and this could be an underestimate (Gaensler 
et al. 2008; Chatterjee et al. 2009; Roberts 2011). 
Therefore we approximate the distance as 2 kpc and 
parameterize it as d = 2d,2 kpc, with a nominal value of 
d,2 = 1.2. In what follows, other fundamental parameters 
for PSR J1816+4510 that are not explicitly cited are 
based on Stovall et al. (2012). 

2. OPTICAL AND ULTRAVIOLET DATA AND ANALYSIS 

2.1. Archival Optical/ UV Data 

We initially identified a potential counterpart to 
PSR J1816+4510 in the USNO-B1.0 survey (Monet 
et al. 2003): the star 1351-0294859 is at 18 h 16 m 35!93, 
+45°10'34!2. This is 0'.'4 from the radio position, reason- 
ably consistent with typical astrometric accuracy from 
the USNO catalog. However, the mean epoch of those 
data is 1974, so a small proper motion could also account 
for some of the difference. The photometry for this source 
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Fig. 1. — Orbital behavior of the radio times-of-arrival (residuals 
at 820 MHz in fjs from Stovall et al. 2012; bottom), Fermi LAT 
0.3-10 GeV count-rate (second from bottom), Swift UVOT u-band 
photometry (relative to the mean; second from top) , and CRTS V- 
band photometry (relative to the mean; top) versus orbital phase, 
repeated twice for clarity. The horizontal lines are for reference 
at the mean magnitude, count-rate, and at residual. The Swift 
data have a scatter consistent with the uncertainties, with x 2 = 
5.4 for 7 degrees-of-freedom relative to a constant model. For the 
Fermi data, x 2 = 20.9 for 15 degrees-of-freedom. The CRTS data 
may show some orbital modulation (x 2 = 29.2 for 14 degrees-of- 
freedom) that shares some of the shape of the Fermi data, but that 
needs to be confirmed. The vertical dashed lines are the times 
of conjunction (eclipse), while the vertical dotted lines show the 
approximate observed limits on the eclipse duration. In the radio 
panel, the thick green segments show the phase region where we 
have observations but the source was not detected. 

is presented in Table 1, where we have assumed uncer- 
tainties of 0.2 mag for the Digitized Sky Survey (DSS) 
photometry. 

The same source is identified in the Catalina Sur- 
veys Data Release 1 (CSDR1; Drake et al. 2009) 
Catalina Real-Time Transient Survey (CRTS). The au- 
tomated software actually identified two sources that 
make up the counterpart: CSS J181635. 9+451033 and 
CSS J181635. 9+451036. These sources are quite close 
(< 2'.'5 apart, which is comparable to the plate-scale 
of the instrument) and they both have the same aver- 
age magnitude of Vcss ~ 18.4 mag, measured with an 
unfiltered detector. There were no images where both 
sources were seen at the same time, and comparing the 
positions of individual detections (rather than the aver- 
age positions in the catalog) it seems they are the same 
source that was split by the photometric pipeline into 
two. The photometry of the combined source (106 mea- 
surements over 6.5 yr, from 2005 May to 2011 Novem- 
ber) is largely consistent with being constant, except for 
8 points. Three of those are clearly when the software 
mis-identified a slightly brighter (Vcss ~ 16.5 mag) star 
about 10" to the South East (visible in Figure 2) as being 
part of this object. The others are not as easy to reject, 
but since they also have Vcss ~ 16.5 mag we think it 
likely that it was another mis-identification or a photo- 
metric artifact; without the images we cannot be certain. 
Excluding those 8 points we have data consistent with a 
constant Vcss = 18.47 mag with root-mean-square vari- 
ations of 0.18 mag. The x 2 relative to a constant model 
is slightly high (156.1 for 97 degrees of freedom), but 
is similar to that for a star of similar brightness 30" 
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Fig. 2.— Optical images of the PSR J1816+4510 field. The ra- 
dio position is indicated with the tick-marks, and the uncertainties 
are dominated by uncertainties in the absolute astrometry of the 
optical data. The large image is the WIYN MiniMo i?-band data, 
showing a 2' X 2' portion. The white band is the gap between the 
MiniMo CCDs. The contours are from the GALEX near-UV im- 
age, and the source GALEX J181635. 9+451034 is consistent with 
a point source. The insets on the right are (top to bottom): the 
WIYN MiniMo R-band data zoomed to show a 30" X 30" por- 
tion and with a gray-scale adjusted to show the faint structure; 
the WIYN MiniMo ij-band data smoothed with a Gaussian kernel 
with cr = 2 pixels (tf'28); and the Swift UVOT M-band image. 

away. There is no evidence for any secular trends in the 
photometry. The uncertainty on the mean magnitude is 
about 2 mmag, but for absolute photometry we transform 
from the unfiltered instrumental magnitude to Cousins 
V by V = Vcss + 0.31(5 - Vf + 0.04 with a scatter of 
0.06 mag 15 . Since this object has very nearly B — V « 
(this assumes that the colors are constant over time and 
orbital phase), we find V = 18.51 ± 0.06, although this 
uncertainty may be somewhat underestimated. 

We then identified the same source in the GALEX 
All-sky Imaging Survey (AIS; Morrissey et al. 2007) 
database. Here, the source is GALEX J181635.9+451034 
and its position is offset by 0'.'4 relative to the radio po- 
sition, which is consistent with the median offset of 1" 
found for GALEX 16 . Again, photometry is presented in 
Table 1. 

In the USNO source catalog, the average density 
of sources brighter than R = 18.4 mag is 6.1 x 
10~ 4 arcsec -2 , so the false association rate given the mea- 
sured source offset is 3 x 10~ 4 , making it very likely that 
we have the correct counterpart. Moreover, the pres- 
ence of a GALEX source at the same position with a 
significant far-UV (FUV; 152 nm) detection makes it es- 
sentially certain: there are only 4 FUV detections with 
magnitudes brighter than 19.9 mag within a 10' radius, 
so the false association rate is 2 x 10~ 6 . As we will see 
below, this source is so bright and blue that it is rather 
unusual, making the chance of finding one within 0'.'5 of 
the radio position by chance extremely low. 

We then identified an observation with the Swift satel- 
lite (Observation ID 00041440003). We see a bright 
source at the radio position in the data from the Ul- 
traviolet and Optical Telescope (UVOT; Roming et al. 
2005). The observation was on 2010 August 04, and 
consisted of 8 separate integrations in the u filter (cen- 
tral wavelength of 3450 A) spread over 11 hours with a 

15 See http://nesssi.cacr.caltech.edu/DataRcleasc/FAQ.html. 

16 See http://galcxgi.gsfc.nasa.gov/docs/galex/Documents/ERO_ 



total integration of 3173 s and 2 x 2 pixel binning. We 
determined both time-resolved and summed photometry 
from these data using Swift data-reduction tools. First, 
we ran the task uvotsource on each separate observa- 
tion (along with respective exposure maps) to measure 
how bright the object was in each individual integration 
(with the 2011 October 31 calibration database). We 
then summed the integrations using uvotimsum and mea- 
sured the summed magnitude using uvotsource, where 
in both cases the source region was a circle with 5" ra- 
dius centered on the radio position and the background 
region is 25" in radius centered near the pulsar but not 
including any visible sources. The best-fit position of the 
source is 18M6'35'/93, +45°10'34"0, or 0'.'12 away from 
the radio position. This is without any additional bore- 
sight correction beyond that applied by the Swift pro- 
cessing. The final detection significance in the summed 
image was 72.8 cr, but our photometry includes the sug- 
gested systematic uncertainty of 0.02 mag in addition to 
the statistical uncertainty. 

In Figure 1 we show the measured Swift and CRTS 
photometry as a function of orbital phase, where the 
Swift observation times have been corrected to the So- 
lar System barycenter using the barycorr task, and the 
CRTS observation times have been corrected to the hc- 
liocenter using the rvcorr task in IRAF (for an 8-hr orbit, 
the differences between helio- and bary-center are negli- 
gible). The Swift data are consistent with no variations, 
with \ 2 = 5.4 for 7 degrees-of-freedom. Each individual 
measurement can largely be considered instantaneous, as 
the exposure times are at most 800 s compared to an ob- 
servation duration of 11 h and an orbital period of 8.7 h. 
The rms scatter of the data about the mean is 0.05 mag. 
The CRTS data have been binned, with between 2 and 16 
individual observations averaged into each point. There 
may be a trend with orbital phase, such that the data 
are slightly (15%) fainter near phases 17 of 0.75, but while 
formally significant (\ 2 — 29.2 for 14 degrees-of-freedom) 
we do not know the level of systematic uncertainties due 
to artifacts and mis-identifications. There were no mea- 
surements in the bin just before the apparent flux mini- 
mum, but without the raw data we cannot say whether 
there were no observations or just no detections. The 
rms scatter of the data about the mean is 0.08 mag. The 
possible trend in the CRTS data is not necessarily seen 
in the Swift data, although it is difficult to be certain. 

2.2. New Optical Data 

We observed PSR J1816+4510 using the Mini-Mosaic 
Imager (MiniMo) on the 3.5-m Wisconsin Indiana Yale 
NOAO (WIYN) telescope (Saha et al. 2000). The 
data were a 300 s exposure in the Harris- R filter on 
2011 March 24 taken shortly before sunrise. Seeing was 
0'.'7, with a plate-scale of 0'.'14 pixel -1 . The data were 
corrected for bias level and flatfielded using standard 
procedures in MIDAS. The image was registered to the 
International Coordinate Reference Frame (ICRF) us- 
ing 130 Two-Micron All-Sky Survey (2MASS; Skrutskie 
et al. 2006) stars, giving fits with rms residuals of 0'.'2 
in each coordinate. We did photometric calibration us- 

17 Note that for pulsars ephemerides, conjunctions occur at 
phases of 0.25 and 0.75, with 0.75 having the pulsar between the 
dataijioMpaptifaima-id-tthc observer. 
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Fig. 3. — Color-magnitude diagram for the sources in the WIYN 
ij-band and Swift u-band images. The source at the position of 
the pulsar is the filled circle. The black line is the approximate 
u = 22.5 detection limit of the UVOT image. We show the results 
of a reddened (Ay = 1 mag) blackbody model for the temperatures 
shown, with a size of O.OO5-R0/1OO pc. A reddening vector showing 
1 magnitude of extinction is also shown. 

ing an observation of the Landolt 98 field (Stetson 2000) 
earlier in the night, determining a zero-point using 22 
stars; we estimate a zero-point uncertainty of 0.03 mag. 
As in the other data, there was a bright source at the 
position of the pulsar, but here the position offset was 
only 0'.'02. We measured the object using sextractor 
(Bertin & Arnouts 1996), with the same settings that 
we used for the standard stars. In addition to the po- 
tential counterpart, we see some faint structure to the 
north-east and south-west, at distances of 2" to 4". This 
may just be some faint stars near the detection limit of 
the image, but they are also somewhat suggestive of an 
Ha bowshock nebula such as that seen around the black- 
widow system PSR B1957+20 (Kulkarni & Hester 1988) 
or around the non-interacting pulsar/ white dwarf binary 
PSR J0437-4715 (Bell et al. 1995). Images of both the 
WIYN and Swift data are shown in Figure 2. 

2.3. Optical/ Ultraviolet Spectral Energy Distribution 

The potential counterpart of PSR J1816+4510 is very 
blue compared to nearby sources. In Figure 3 we show a 
color-magnitude diagram using the i?-band and w-band 
data, where the counterpart is roughly 1 magnitude bluer 
than the field sources. Just the u and R data indicate a 
rather hot blackbody, although since the reddening vec- 
tor is roughly parallel to the track of a blackbody the 
temperature and reddening are degenerate. However, 
some temperature above 10,000 K is required. 

We fit all of the optical/ UV photometry to determine 
the spectral energy distribution (SED) of this source. For 
the data in the Vega system, we used zero-point fluxes 
from Bessell, Castelli, & Plez (1998). We then convolved 
various model SEDs with filter transmission curves and 
compared the resulting fluxes with those derived from the 
data. For the DSS and CSS data we assumed standard 
Johnson BRI filters, and this is clearly a simplification, 
but the large uncertainties for DSS mean that those data 
have modest weight. For the WIYN data we used a fil- 
ter curve from Kotulla et al. (2009); for the Swift data 
we used a response file from the Swift web site 18 ; the 



GALEX filters were from the COSMOS web site 19 . We 
used a nominal extinction curve from Cardelli, Clayton, 
& Mathis (1989) and O'Donnell (1994), with a reddening 
ratio Rv = 3.1. 

Our first model was a reddened blackbody. We got a 
good fit, with x 2 = 6.4 for 5 degrees of freedom (and 
this includes the poorly calibrated DSS data). The best- 
fit model had T cS = 18,000 K, A v = 0.77 mag, and 
size R = O.li?© at a nominal distance of 2kpc, but a 
wide range of solutions had similarly good fits (Figure 4) 
with larger temperature requiring larger extinctions and 
smaller radii; as in Figure 3 this is largely the result of 
the blackbody model between the u and R data being 
parallel to the reddening vector. The size is largely de- 
termined from the WIYN observation, with 
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along the best-fit locus. 

We then fit model stellar atmospheres from Kurucz 
(1993). We used models with gravity 10 4 cms~ 2 , al- 
though the results were not sensitive to this. The best- 
fit model was slightly hotter than the best-fit blackbody 
(21,000 K, with A v = 0.92 mag and R = 0.12i? Q at a 
distance of 2kpc), and had a slightly worse \ 2 (H-5), 
but uncertainties in the extinction law below 250 nm 
(Cardelli et al. 1989) could change the result; the most 
discrepant point was the GALEX FUV observation. The 
best-fit region in the (Av,T e g) plane is very similar 
to that of the blackbody (Figure 4). We also tried a 
15,000 K white dwarf atmosphere model (a DA white 
dwarf with hydrogen on the surface, which is expected 
for such a hot star; Hansen & Liebert 2003), kindly sup- 
plied by D. Koester. The fit was similar to that of a 
15,000 K main-sequence star since the main differences 
(the stronger Balmer absorption sequence in the white 
dwarf) are not easily distinguished with the available 
DSS [/-band photometry; the white dwarf model also 
had trouble fitting the GALEX FUV point. 

An absorbed power-law spectrum such as that of an ac- 
tive galactic nucleus does not fit (% 2 = 34.7 for 5 degrees 
of freedom). While the spin-down and radio emission of 
PSR J1816+4510 indicate that no accretion disk is likely 
present (cf. Archibald et al. 2009), given the unusual na- 
ture of the optical/UV emission it may be worth consid- 
ering whether the high temperature that we measure is 
from an accretion disk around the pulsar, where ultra- 
violet emission is common (although this would tend to 
also have bright X-ray emission which we do not see, the 
X-rays could be variable). The SED that we measure 
is consistent with a single-temperature blackbody, not 
a multi-temperature model as is usually used to model 
accretion disks (Shakura & Sunyaev 1973; Vrtilek et al. 
1990). Therefore, while we cannot rule out such a model 
using photometry alone, we consider it unlikely Spec- 
troscopy should be definitive. 

We note that the faint, diffuse emission seen in the R- 
band image could have contaminated some of the lower- 



http:/ /heasarc. nasa.gov/docs/swift/proposals/swift_responses. html. 19 http:/ /www. astro. caltech.edu/-capak/cosmos/filters/. 
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Fig. 4. — The results of fitting to the broadband photometry from 
Table 1. We plot the best-fit SEDs, with the blackbody model 
(black) and model atmosphere (red). The solid curves have had 
extinction applied, while the dashed curves are without extinction. 
The black squares and red diamonds are the respective model SEDs 
convolved with the filter passbands, while the circles are the data. 
The inset shows contours of x 2 vs - extinction Ay and effective 
temperature T ef j for the blackbody models (black) and model at- 
mospheres (red). Contours are at 1-, 2-, and 3-<r confidence on 
the joint fit. The stars represent the best-fit parameters for each 
model. For comparison we also plot the best-fit blackbody with 
Ay = 0.2 mag (green). 

resolution data (in particular Swift and GALEX). How- 
ever, at least at i?-band the brighter spot (to the south- 
west) is about 100 times fainter than the star, so as- 
suming a similar spectrum (which is conservative for the 
bluer bands, as the emission is probably either Balmer 
dominated or stellar) it is unlikely to be important at 
the > 1% level. We tried subtracting the stellar point- 
spread-function (PSF) at the position of the counterpart, 
but do not see any significant residual emission beyond 
that visible in Figure 2. 

3. X-RAY AND 7-RAY DATA AND ANALYSIS 
3.1. Swift X-ray Data 

In the 2.8 ks photon-counting observation with X-ray 
Telescope (XRT; Burrows et al. 2005) counts were de- 
tected in a circle with radius 20" around the radio po- 
sition. Millisecond pulsars typically have a combination 
of X-rays from thermal (from hot polar caps) and non- 
thermal (either magnetospheric, or from a shocked pul- 
sar wind) spectra (e.g. Zavlin 2007). The thermal com- 
ponents have blackbody temperatures with kT w 0.1- 
0.2 keV and luminosities of - 10" 3 £ (Zavlin 2007; Pos- 
senti et al. 2002). For more energetic pulsars the non- 
thermal components dominate, and these are typically 
fit as power-laws with photon indices T « 1.5 and 
again luminosities of ~ 10~ 3 i?. Based on a power- 
law with photon index of 1.5, we set a 2 a limit of 
< 3 x 10~ 15 ergs -1 cm -2 for the unabsorbed 0.5-10 keV 
flux from PSR J1816+4510 for absorption column den- 
sities iVn in the range of 10 20 cm _2 -10 21 cm~ 2 (cor- 
responding to Ay = 0.06 — 0.6 mag, based on Pre- 
dehl & Schmitt 1995). For a measured spin-down lu- 
minosity of E = 5 x 10 34 ergs _1 (although see below 
for possible corrections to this), our limit then corre- 
sponds to ix,non-th ^ 3 x lO -4 ^ 2 ^. We can do the 
same computation for a thermal spectrum, with un- 
absorbed flux limits of (9 — 17) x 10 -15 ergs -1 cm~ 2 
(0.5-2 keV) for a blackbody with kT = 0.15 keV and 



cm 

< 2 x lQ- 4 d^' z E. Both of these 
efficiencies are low but not outside the observed range 
(Zavlin 2007; Ransom et al. 2011), and suggest that the 
X-ray flux may only be slightly below the Swift limit. 

3.2. Fermi Data 

The radio position of PSR J1816+4510 matches almost 
exactly with a source from the Fermi Large Area Tele- 
scope Second Source Catalog (2FGL; The Fermi-LAT 
Collaboration 2011). The source 2FGL J1816.5+4511 
(1FGL J1816.7+4509 from the first year catalog) is 1.4' 
away from PSR J1816+4510, with a position uncertainty 
of w 5' in radius. It is listed as having a power-law spec- 
trum with photon index T = 2.11 ± 0.08 (N E oc £~ r ), 
and 0.1-100 GeV flux of (15.3±1.8) x 10~ 12 ergs -1 cm -2 . 

However, there is another year of data available in the 
Fermi archive, and we wished to do further spectral anal- 
ysis and look for pulsations. Therefore we analyzed data 
from the Fermi Large Area Telescope (LAT; Atwood 
et al. 2009), including events from 2008 August 05 to 
2012 January 19. We followed standard procedures 20 in 
filtering events, selecting those with event class 2 within 
a 10° radius around PSR J1816+4510, with energies be- 
tween 0.2-10 GeV (to avoid the poor point-spread func- 
tion at the lowest energies; the pulsar did not appear to 
be detected there anyway) and zenith angles < 105°. 
We computed the spectrum using an unbinned likeli- 
hood analysis 21 , including the contributions of sources 
out to a radius of 17° from the 2FGL catalog as well 
as isotropic and Galactic background models appropri- 
ate for pass 7 data (models iso_p7v6source.txt and 
gal_2yearp7v6_v0.f its) with the P7S0URCE_V6 instru- 
ment response, although we held most of the source pa- 
rameters fixed at their catalog values with the excep- 
tions of those sources within 8° of PSR J1816+4510 
and the diffuse background normalizations. Photons 
are significantly detected between 0.5 GeV and 5 GeV. 
For 2FGL J1816.5+4511 we find a good fit with a 
power-law model with photon index T — 2.20 ± 0.07 
and normalization (1.5 ±0.1) x 10~ 12 cm -2 s -1 MeV -1 
at 1.15 GeV giving an integrated 0.1-100 GeV flux of 
(19.6 ± 1.5) x 10~ 12 ergs -1 cm" 2 , with a Test Statistic 
of 404 (i.e., roughly a 20 a detection); this is similar to 
the result from the 2FGL catalog. We do not incorpo- 
rate any systematic uncertainties related to calibration 
errors. 

We repeated the spectral fit with a power-law modi- 
fied by an exponential cutoff, N E oc E~ v e~ E l Ea . The 
meager energy range with significant detections meant 
that the cutoff could not be strongly constrained, but 
wc find T = 2.0 ± 0.1, normalization (7.9 ± 0.8) x 
at 0.55 GeV, and cutoff energy 
7.5 ± 4.0 GeV; the integrated flux was (15 ± 3) x 
Formally this fit was statistically 
equivalent to the pure power-law fit, and other local min- 
ima were also possible depending on where the fit was 
started. We note that the parameters we find are out- 
side the range of most millisecond pulsars, with F and 
E c both higher than are typical. Much of this comes 

20 See http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/data_preparation 

21 See http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/likclihood_tutorii 
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from the highest energies we included in our fit (4-6 GeV; 
Fig. 5). Without this bin, a fit with more typical values 
(r w 1.5, E c Ri 3 GeV) is acceptable. We show the fits in 
Figure 5, where 7-ray fluxes were determined from mod- 
eling the flux in each energy bin as a single power-law 22 
using the contributed task likeSED. We therefore urge 
caution in interpreting the 7-ray spectrum. 

For the pulsation search, after the initial event filter- 
ing, we assigned phases to all of the events using the best- 
fit radio ephemeris using the Fermi plugin for tempo2 23 
(Hobbs, Edwards, & Manchester 2006). We detected 
pulsations using the initial radio ephemeris, but given 
the longer time span of the Fermi data compared to the 
radio data (1262 days vs. 320 days) we were able to re- 
fine the radio ephemeris (in particular the spin-down), 
as discussed in Stovall et al. We then used the refined 
ephemeris to update the event phases. Selecting the 632 
events < 0.65° from the radio position and with energies 
between 0.3 GeV and lOGev (optimizing those parame- 
ters for the pulse amplitude, as in Ransom et al. 2011), 
we see very significant pulsations, with a H-test statistic 
(de Jager, Raubenheimer, & Swanepoel 1989) of 64.4 for 
11 harmonics (false alarm probability of 4 x 10 -23 ). With 
this solution we see clear pulsations in the binned light 
curve with x 2 = 157 for 19 degrees-of-freedom. The pul- 
sations have two sharp peaks separated by slightly less 
than one half of a cycle (Fig. 6) similar to the radio pulse. 
Selecting events from phases 0.32-0.52 (where the pulse 
is at a minimum), we see a radial profile that is consis- 
tent with being flat in terms of counts per unit area out 
to 2°, with an implied background rate over all phases 
of 237 ± 38deg -2 (the horizontal line in Fig. 6). This 
would mean that our lightcurve is consistent with being 
100% pulsed. Averaging over pulse phase, the count- 
rate was roughly constant as a function of orbital phase 
(Fig. 1), with x 2 — 20.9 for 15 degrees-of-freedom. Over 
the two years of the 2FGL catalog the lightcurve was like- 
wise consistent with being constant on month timescales 
(x 2 = 22.1 for 23 degrees-of-freedom, based on the 2FGL 
variability index). 

4. DISCUSSION 

Based on the Fermi flux of 20 x 10~ 12 erg cm -2 s _1 
(0.1-100 GeV) we find a 7-ray luminosity of 1 x 
lO 34 ^ 2 , ergs" 1 . We determine the 7-ray efficiency by 
comparing this with E and find ?7 7 = L 7 /E = 0.19^ as- 
suming no geometric beaming correction; this is compa- 
rable with that found for purely magnctospheric emission 
from millisecond pulsars with Fermi (Abdo et al. 2009; 
Ransom et al. 2011). We can set a weak upper limit to 
the distance by requiring L 1 < E, which gives d 2 < 2.3 
(d < 4.6 kpc). Given the highly pulsed 7-ray lightcurve, 
it is possible that all of the emission could be magne- 
tospheric in origin like for most millisecond pulsars, but 
some might still be related to intra-binary shocks such as 
in the PSR B1259-63 system (Abdo et al. 2011, which 
is not a millisecond pulsar but has a hot companion like 
that seen here) , especially if there are contributions from 
inverse Compton scattering off the UV photons present 
in both cases. 

22 See http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/python. 

23 See http:/ /fermi. gsfc.nasa.gov/ssc/data/analysis/scitools/pulsar. 
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Fig. 5.— The SED of PSR J1816+4510 and its presumed com- 
panion, from the optical to 7-rays. We show the optical/UV pho- 
tometry along with the best-fit blackbody model, both corrected 
for reddening with Ay = 0.2 mag. The X-ray upper limit is based 
on the Swift XRT non-detection. The 7-ray points were derived 
from fitting a power-law model to each individual energy bin, and 
we show the best-fit single power-law and cutoff power-law models 
along with ±lcr uncertainties. 

In Figure 5 we plot the spectral energy distribution 
from optical to 7-rays. Energetically, the optical/UV are 
almost as important as the 7-rays, which would make 
it difficult for them to both ultimately come from E 
(as the 7-rays already require a substantial fraction of 
E), supporting a hot companion which radiates on its 
own. Changing the extinction to a lower value such as 
Ay =0.2 mag (see below) reduces the total optical lumi- 
nosity somewhat, but it is still substantial. The simple 
power-law fit to the 7-rays exceeds the X-ray upper limit 
(much as in Durant et al. 2011), but the cutoff power- 
law does as well, so with either model we might need 
a spectral break between 1 keV and 100 MeV. However, 
our spectral fitting only had a limited number of counts, 
and we did not include systematic uncertainties related 
to instrumental calibration. The apparent discrepancy 
between the 7-ray emission and the X-ray upper limit 
may only be a consequence of the spectral fit; the ra- 
tio of ~ 10 3 between the 7-ray and X-ray luminosities is 
reasonable for other millisecond pulsars (Ransom et al. 
2011; Takata, Cheng, & Taam 2012). 

In what follows, we primarily assume that the opti- 
cal/UV emission come from a single photosphere that is 
the companion of the pulsar in a binary system. How- 
ever, without a spectrum we cannot exclude contribu- 
tions from shocked plasma — this might help explain 
some of the slightly discrepant UV data-points or the 
somewhat high extinction (see below). 

Given our extremely likely detection of the optical/UV 
companion of PSR J1816+4510, and the identification of 
the pulsar at 7-ray energies, we consider how this sys- 
tem fits among the known recycled pulsars with low-mass 
companions. Some of the implications of this system 
are common to a wide range of similar systems, but the 
uniquely hot temperature may point to cither a fortu- 
itous detection of a short-lived evolutionary state or a 
different evolutionary path. 
At the nominal distance, our optical photometry im- 
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Fig. 6.— The pulsed 7-ray light curve of 2FGL J1816. 5+4511, 
determined from all LAT data in the 0.3— lOGeV range up to 
2012 January and repeated twice for clarity. We used the ra- 
dio ephemeris to initially identify the pulsar, but then refined the 
ephemeris using the Fermi data since they have a longer time span. 
The dashed horizontal line is the approximate background level 
based on the phases indicated by the vertical dotted lines. The red 
trace is the radio pulse profile, based on 820-MHz data from the 
Green Bank Telescope (Stovall et al. 2012), arbitrarily scaled. 

ray properties but is unlikely (it would require L 1 ss AE, 
but this constraint is based on isotropic emission). If the 
optical source is an unrelated object, it would either be 
a nearby white dwarf at < 300 pc (Fontaine, Brassard, & 
Bergeron 2001 24 ; we take it as a O.6M carbon/oxygen 
white dwarf) or a main sequence star (B3-5) at a dis- 
tance of 70-100 kpc (Cox 2000); while the former is pos- 
sible, main sequence stars with M > 3Mq would not be 
expected at such distances. 

Our best-fit value for the extinction Ay is 0.5-1.0 mag. 
This is larger than the largest value in this direction 
(ss 0.2 mag) determined by Drimmel, Cabrera-Lavers, 
& Lopez-Corredoira (2003), and Schlegel, Finkbeiner, 
& Davis (1998) give a similar result. Given the limi- 
tations of our fit, values of < 0.2 mag are not excluded 
but would imply effective temperatures of ss 12,000 K 
(Ay = 0.2 mag increases \ 2 by 2.2 over the best-fit 
value for the blackbody fit), and we plot the best-fit 
Ay — 0.2 mag blackbody for comparison in Figure 4. 
To evaluate the likely extinction, we determined our own 
run of extinction with distance by examining all of the 
2MASS stars within 1° and finding the red clump (Drim- 
mel et al. 2003; Durant & van Kerkwijk 2006). While 
we cannot determine the extinction as close as 2 kpc 
(there are not enough stars), we measure extinctions of 
ss 0.6 mag for distances > 5 kpc, which is reasonably 
consistent with our SED fitting. At this Galactic lati- 
tude (b = +24.7°) much of the extinction will be close 
to the Sun, so the value measured at 5 kpc should be 
applicable to closer objects. We note that our value 
of Ay > 0.5 mag is actually consistent with the mea- 
sured DM, assuming an ionized fraction of 10% and the 
usual (Predehl & Schmitt 1995) conversion between ex- 
tinction and hydrogen column density, but the hydrogen 
column density interpolated from H I maps 25 is lower, 
4 x 10 20 cm~ 2 . Optical and X-ray spectroscopy can hope- 
fully narrow down the possible range of the extinction. 

24 Also see http:/ /www. astro. umontreal.ca/~bergeron/CoolingM< 

25 Using http://asc.harvard.edu/toolkit/colden.jsp. 



With an minimum eclipse duration of 7% of the or- 
bit at 820 MHz (Stovall et al. 2012), the eclipsing ra- 
dius is w 0.5-R©. This is similar to the Roche lobe 
radius, suggesting that some of the eclipsing material 
may be gravitationally bound to the companion star, al- 
though the tail of delayed times-of-arrival (TOAs) ex- 
tends to larger radii. The maximum delay observed 
at 820 MHz (about 400 /is) implies an excess dispersion 
measure of 0.06 pc cm -3 , or an electron column den- 
sity of N e ss 2 x 10 17 cm~ 2 . If this material is dis- 
tributed over 0.5-R©, we would have an electron density 
n e ps 6 x 10 6 cm~ 3 . Assuming that the material is mov- 
ing at the escape velocity, we estimate a mass-loss rate 
of M ~ 10~ 13 M Q yr~ 1 , so the companion would not be 
substantially diminished over a Hubble time (similar to 
Stappers et al. 1996). Such mass-loss rate is actually 
comparable with expectations for radiative winds from 
more massive sdB stars with (presumably) similar sur- 
face gravities 26 (Vink & Cassisi 2002; Unglaub 2008). If 
this and not ablation (which would only require 0.1% of 

E) is the origin of the mass loss, then the low gravity and 
high temperature of the companion seem to be necessary 
components for the ionized gas eclipses as winds cease 
for gravities > 10 6 cms -2 and temperatures < 20,000 K 
(Unglaub 2008). While minor in terms of mass loss the 
winds might be nonetheless important evolutionarily in 
altering the diffusive equilibrium (Unglaub & Bues 1998) 
and hence the atmospheric appearance and onset of shell 
burning. 

The measured effective temperature of > 15,000 K is 
far hotter than the companion to any known black- widow 
or redback system (typically < 6000 K; van Kerkwijk 
et al. 2005; Pallanca et al. 2010; Breton et al. 2012, in 
prep.; C. Bassa 2011, pers. comm.) by a factor of al- 
most 3. Pulsars with hot white-dwarf companions are 
known (e.g., 8,550 K for PSR J1012+5307, 10,090 K for 
PSR J1911-5958A, and 15,000 K for PSR B0820+02; 
van Kerkwijk et al. 2005 and references therein), but they 
do not have broad eclipses like those we see here. We can 
then address the peculiarity of this system in two ways: 
(1) if this is an interacting binary system, what would be 
the consequences of it being so hot, and (2) why is it so 
hot. 

For the first question, if we treat PSR J1816+4510 as 
an interacting binary regardless of origin, it is not sur- 
prising that we do not see any modulation of the it-band 
lightcurve: the additional energy deposited by the pul- 
sar has an equilibrium temperature < (£ , /47ra 2 cr) 1 / 4 w 
7000 K (for an efficiency of 100%, while typical efficien- 
cies are closer to 10%; Breton et al. 2012, in prep.). 
Note that this ascribes all of the observed spindown 
P = 4.1 x 10 _20 ss _1 to magnetic dipole radiation and 
not to secular (i.e., Shklovskii 1970) acceleration; with a 
typical millisecond pulsar velocity of v — 100viookms~ 
the secular P would be 1.7 x lO -21 ^^^ 1 s s _1 or 4% 
of the measured value, so it is likely that the E value 
we use above is close to correct (corrections due to dif- 
ferential Galactic acceleration are even smaller; Nice & 
Taylor 1995). Since the implied temperature is so low 

26 Although we are outside their nominal luminosity range and 
near the low end of the temperature range usually considered, the 
surface gravities are similar so the winds are likely to be similar. 
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compared to the observed temperature, we would expect 
the illuminated side to be < 5% brighter than the dark 
side (for an effective temperature of 15,000 K), which is 
the level of the observed scatter in the Swift photome- 
try; if the orbital modulation of the CRTS data is real, 
then it is difficult to understand its amplitude. We note 
that our SED fitting assumes that the data arc con- 
stant with time and are defined by only a single model, 
and these are not necessarily valid assumptions. Even 
aside from the lack of strong orbital modulation, there 
could also be secular/state changes such as those seen in 
PSR J1023+0038 (Archibald et al. 2009). However, the 
limit on secular evolution from the CSDR1 data (which 
span the times of the GALEX, Swift, and W1YN Pho- 
tometry) suggests that most of the data are consistent 
with a single model. In the future, single-epoch pho- 
tometry should be able to resolve this question. The 
nominal radius of the companion is well within the typ- 
ical range for black- widow/redback companions in the 
field (sources in globular clusters often are considerably 
> Rq] Ferraro et al. 2001; Cocozza et al. 2008; Pallanca 
ct al. 2010), which may be some clue to the formation 
mechanism. However, unlike some of those systems, it is 
likely a factor of several smaller than the Roche lobe, so 
the companion may not be significantly distorted. 

It may be that the O.5i?0 eclipse duration is set not by 
the Roche lobe radius but by the intrabinary shock radius 
between the £?-driven wind and the companion (Arons 
& Tavani 1993), where the shock represents the equi- 
librium between the relativistic wind of the pulsar and 
the presumed radiation-driven wind from the companion. 
Assuming an electron density of ~ 10 7 cm~ 3 and pure 
hydrogen composition at a radius of i? s hock ~ 0.5i? Q , 
the ram pressure pv 2 /2 is 0.03 dyne cm -2 for a wind at 
the escape velocity, which is a factor of ~ 200 less than 
E I \4irc(a — -R s hock) 2 ); if this model is valid, then we 
must consider that the wind might be moving faster than 
the escape speed (Phinney et al. 1988; Unglaub 2008), 
our density might be too low (in particular the mate- 
rial might be clumpy), or possible only a fraction of E 
participates in the shock (Stappers et al. 2003). 

For the second question, we can ask why the com- 
panion would be so hot. First, it is possible that 
PSR J1816+4510 has a normal cool companion, but that 
the optical/UV flux comes from another source. This 
could be a star, either as part of a triple system or an un- 
related object. Having an unrelated object seems highly 
unlikely given the positional coincidence, and an asso- 
ciation between the radio, 7-ray, and optical sources is 
the most likely explanation, but without optical modula- 
tion to confirm we cannot be certain. A triple system can 
largely be ruled out by the 3-yr span of the Fermi timing, 
as those data are consistent with only the 8-hr binary. As 
mentioned before, emission from shocked plasma is also 
possible. 

Without knowing the surface gravity (and hence hav- 
ing some idea of how degenerate the companion is) , any 
inferences about why the companion is so hot are dif- 
ficult. For helium-core white dwarfs in the mass range 
considered here (< 0.2M Q ) burning of a thick shell of 
hydrogen (Webbink 1975) can keep the sources hotter 
for considerably longer (Gyr; Sarna et al. 2000; Panei 
et al. 2007; Steinfadt, Bildsten, & Arras 2010) than stan- 



dard cooling would allow (cf. Lorimer et al. 1995), but 
the temperatures tend to be < 10,000 K; more massive 
white dwarfs can stay at > 10 4 K for longer, and un- 
stable hydrogen flashes can also push the temperature 
above 10 4 K for more massive white dwarfs temporar- 
ily. We note, though, that the recently discovered bi- 
nary SDSS J065133.33+284423.3 (Brown et al. 2011) has 
a helium-core white dwarf at a similar temperature to 
what we find, and this source is somewhat more mas- 
sive (0.25M Q ) than the expected burning limit, although 
the limit is metallicity-dependent and the contribution of 
tidal heating in this system is unknown. If the compan- 
ion of PSR J1816+4510 were a normal white dwarf we 
would expect PSR J1816+4510 to be on the low side of 
the possible distances (radii of < 0.05i? Q are typical at 
these masses; Panei et al. 2007) with an age of < 0.5 Gyr. 
In this scenario, PSR J1816+4510 would be the first pul- 
sar/white dwarf system with ionized-gas eclipses. How- 
ever, it could also be that the companion is still hotter 
and younger still if it has not yet equilibrated from a 
common-envelope evolutionary phase (Paczynski 1976) 
or Roche-lobe overflow that stripped away the outer lay- 
ers of the companion, leaving it hot, large, and in a 
circular orbit (e.g., Driebe et al. 1998). For low- mass 
(< 0.2 Mq) white dwarfs, the residual hydrogen burn- 
ing typically results in luminosities of < O.IL© (Panei 
et al. 2007; Steinfadt et al. 2010), and therefore it might 
be that the luminosity we see is dominated by gravita- 
tional contraction (but see Driebe et al. 1998). If that 
is true, we estimate a thermal timescale of 2<i^~ 3 Myr 
for an effective temperature of 18,000 K and consider- 
ing the whole star — the timescale would be 2d,2 3 kyr if 
we only are concerned with a typical envelope of 1O _3 M 
(Panei et al. 2007) — so if the companion is larger than a 
typical low-mass white dwarf then it is likely extremely 
young, and we might be seeing a newly-born millisec- 
ond pulsar slightly after the transitional phase found by 
Archibald et al. (2009). The hot/large state could also 
be a result of a recent hydrogen flash. This source then 
resembles the rather young WD/sdB stars HD 188112 
(T cB =21,500K, M = 0.24M Q , R = 0.1 i? ; Hebcr 
et al. 2003) or GALEX J1717+6757 (T cff =14,900 K, 
M = 0.19 M Q , R = 0.1 R Q ; Vennes et al. 2011) which 
are thought to be progenitors of more typical helium-core 
white dwarfs. Both of these scenarios (mass stripping or 
hydrogen flash) are possibilities for the hot, bloated white 
dwarfs seen by Kepler (van Kerkwijk et al. 2010; Carter, 
Rappaport, & Fabrycky 2011; Breton et al. 2012), which 
the companion to PSR J1816+4510 also resembles, al- 
though the presence of a neutron star instead of a main 
sequence primary would require a different evolution. 
Again, spectroscopy to determine surface gravity and el- 
emental abundances should be definitive. Finally, it is 
possible that, despite the eclipses the system is closer 
to face-on than edge on, as for inclinations < 30° the 
companion mass is > 0.5M Q like other hot white dwarf 
companions. This would be contrary to the theoreti- 
cal companion mass vs. orbital period relation (Tauris & 
Savonije 1999) and would make the radius even stranger, 
but there is at least one known outlier from the compan- 
ion mass vs. orbital period relation with a substantially 
higher mass than expected (Hessels et al. 2005; Demorest 
et al. 2010). 
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The bright counterpart makes optical astrometry 
within the reach of ground-based telescopes, at least for 
determining a proper motion. While we see no defini- 
tive proper motion comparing the Swift and WIYN data 
(taken 0.6 yr apart), individual ground-based images can 
determine the relative position of the pulsar to < 10 mas 
in only a few minutes. We expect a proper motion 
fj, = 10-uiooc^ 1 masyr -1 , so assuming adequate calibra- 
tion this can be measured in a year or two. We could 
then compare this against any radio proper motion de- 
termined from timing, which would further establish 
whether or not the optical source is indeed the compan- 
ion of the pulsar. 

5. CONCLUSIONS 

We have discovered the almost certain optical/UV 
counterpart of the newly-discovered millisecond radio/7- 
ray pulsar PSR J1816+4510. The radio/7-ray properties 
of PSR J1816+4510 appear much like most energetic 
eclipsing pulsars discovered recently (Roberts 2011), 
which supports the use of Fermi and low-frequency ra- 
dio observations to find energetic recycled pulsars. Some 
aspects of PSR J1816+4510's radio properties appear 
unique, in that the plausible size of the eclipsing region 
seems to be contained in its companion's Roche lobe, but 
such inferences depend on the inclination as well as the 
observing frequency of the radio data. 

While we still cannot constrain all of its parameters 
uniquely, the optical/UV properties of PSR J1816+4510 
appear more robustly unique, with an effective tempera- 
ture likely at least 3 times higher than any other black- 
widow or redback system. In fact, the companion may be 
the brightest low- mass (i.e., not a B star) optical 27 com- 
panion to any pulsar (van Kerkwijk et al. 2005). This 
is largely because of the high temperature rather than 
a small distance or large size. The high temperature 
presents a number of puzzles and opportunities. Depend- 
ing on the radius, it may be that the companion is rather 
young, and that we are seeing the sources only shortly 
after its envelope was stripped away. Phase-resolved pho- 
tometry and spectroscopy will be important to determine 
the orientation and mass function of the system (van 
Kerkwijk et al. 2011), and radio astrometry can help 
constrain the radius of the companion. Long-term op- 
tical monitoring may be able to detect cooling after a 
recent stripping or burning episode. Given how bright it 



is, modulation of the companion may be detectable from 
a number of mechanisms. Orbital motion may be visi- 
ble through Dopplcr boosting (van Kerkwijk et al. 2010; 
Shporer et al. 2010), which is expected to produce mod- 
ulation of ±0.3%, while ellipsoidal modulations (±0.1%) 
could help constrain the mass ratio and/or the degree of 
Roche-lobe filling (van Kerkwijk et al. 2010; Carter et al. 
2011; Breton et al. 2012). At the same time, the high 
temperature and small size compared to the Roche lobe 
mean that the photocenter may be much closer to the ge- 
ometric center of the star (cf. van Kerkwijk et al. 2011) 
facilitating (along with its brightness) measurement and 
interpretation of the radial velocity curve, and with it 
enabling an accurate mass measurement for this neutron 
star. 
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